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Abstract. The phenomenon of partial equivalence of statistical ensembles is illus- 
trated by discussing two examples, the mean-field XY and the mean-field spherical 
model. The configurational parts of these systems exhibit partial equivalence of 
the microcanonical and the canonical ensemble. Furthermore, the configurational 
microcanonical entropy is a smooth function, whereas a nonanalytic point of the 
configurational free energy indicates the presence of a phase transition in the canonical 
ensemble. In the presence of a standard kinetic energy contribution, partial equivalence 
is removed and a nonanalyticity arises also microcanonically. Hence in contrast to the 
common belief, kinetic energy, even though a quadratic form in the momenta, has a 
non-trivial effect on the thermodynamic behaviour. As a by-product we present the 
microcanonical solution of the mean-field spherical model with kinetic energy for finite 
and infinite system sizes. 



PACS numbers: 05.20.Gg, 05.20.-y, 05.70.Fh 
1. Introduction 

Different statistical ensembles, like the microcanonical or the canonical ones, describe 
different physical situations. Under certain conditions, statistical averages calculated 
in different ensembles are the same in the thermodynamic limit, and this situation is 
referred to as ensemble equivalence. Contrary to the widespread belief among physicists, 
ensemble equivalence does not hold in general. Prominent counterexamples are systems 
with long-range forces, like gravitational or unscreened electrostatic interactions. 

An in some sense intermediate case between equivalence and nonequivalence is 
the only recently described partial equivalence of ensembles. When two ensembles are 
properly nonequivalent, there is a set of values of the control parameter in one ensemble 
that does not correspond to any value of the control parameter in the other ensemble. 
To give an example, there may be forbidden energy values, i. e., energy values for which 
a microcanonical description is perfectly legitimate, but which cannot be realized as 
canonical equilibrium average energies for any value of the temperature. In a situation 
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of partial equivalence, a whole set of values of the control parameter in one ensemble 
corresponds to a single value of the parameter in the other ensemble: the same average 
energy value of the energy is obtained for a whole set of temperatures, or vice versa. 

In the present paper, we focus on a particular kind of partial equivalence between 
the canonical and the microcanonical ensemble, where an interval of temperatures 
corresponds to a single value of the energy. Only a few examples of this kind of partial 
equivalence have been discussed in the literature so far. A crucial ingredient of all 
examples known is a bounded above Hamiltonian function, and for this reason the 
absence of a standard kinetic energy term, i.e., of a quadratic form in the momenta, is 
mandatory. The main goal of the present article is to discuss the phenomenon of partial 
equivalence as well as the way full equivalence of ensembles is restored by adding a kinetic 
energy term to the Hamiltonian function of the system. We observe that adding such a 
standard kinetic energy may have a nontrivial effect on the thermodynamic behaviour 
of the system. 

The paper is organized as follows. After recalling in section [2] the definitions 
of nonequivalence, partial equivalence, and full equivalence of ensembles, we discuss 
two examples, the mean-field XY and the mean-field spherical model, in sections [3] 
and HI Considering the potential energy of these models alone, both systems exhibit 
partial equivalence of the microcanonical and the canonical ensembles. Furthermore, 
the configurational microcanonical entropy is a smooth function, whereas a nonanalytic 
point of the configurational free energy indicates the presence of a phase transition in 
the canonical ensemble. Adding a standard kinetic energy term to the Hamiltonian, 
as shown in section HJ equivalence of ensembles is observed and phase transitions are 
indicated by the presence of nonanalytic points in both, the microcanonical entropy 
and the canonical free energy. This result is at variance with the common belief that 
a standard kinetic energy yields only a trivial term in the thermodynamic potentials. 
We conclude with a discussion of the results and their implications in section [5j In 



Appendix A we supplement these results by pointing out that not only the addition of 



a kinetic energy term, but also a coupling to an external magnetic field may remove 



partial equivalence of the microcanonical and the canonical ensemble. In Appendix B 



we report the exact microcanonical solution, at finite and infinite N, of the mean-field 
spherical model with a kinetic energy term. 



2. (Non) equivalence and partial equivalence of statistical ensembles 

The various statistical ensembles model different physical situations: The microcan- 
onical ensemble describes an energetically isolated system with a given number of 
particles, confined to a fixed volume. The canonical ensemble models a system in 
thermal equilibrium, i. e., exchanging energy with the environment. The grand-canonical 
ensemble is the appropriate choice when a system is in thermal equilibrium and the 
number of particles is not fixed. These are the ensembles most widely used, but 
one may also consider further statistical ensembles, like those where pressure is held 
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constant instead of volume, and many more. In general, the ensemble averages of a 
given observable yield different values, depending on which ensemble is used. However, 
in the thermodynamic limit (in case it exists), averages from different ensembles may 
coincide, and in this case one speaks of ensemble equivalence. Under suitable conditions 
on the interparticle interactions, ensemble equivalence is guaranteed to hold pQ. In 
physical terms, these conditions are satisfied by short-ranged, bounded below potentials 
with a hard- or soft-core short- distance repulsion describing intermolecular forces, like 
the Lennard- Jones potential. In condensed matter systems (which are the common 
objective of statistical physics) these conditions are typically satisfied, and for this reason 
ensemble equivalence is taken for granted in many treatises on statistical physics. 

In recent years, however, nonequivalence of ensembles has attracted a significant 
amount of attention. What has been known among astronomers for almost a century, 
i.e., that long-range forces like gravitational or unscreened electrostatic interactions 
yield different statistical averages when treated canonically or microcanonically, has 
been realized only recently within the statistical mechanics community. For such long- 
range systems, however, the energy is a nonextensive quantity, and a thermodynamic 
limit in the usual sense does not exist. More recently it became clear that also in systems 
where the energy is extensive and a proper thermodynamic limit exists, one can find 
ensemble nonequivalence. The main property of these systems is that the energy, though 
extensive, is nonadditive [2], usually due to the long-range nature of the interactions. 
The simplest examples of systems showing nonequivalence of ensembles are spin models 
with long-range interactions [31 HI El El IE] , other examples include self-gravitating 
systems [91 [101 E] an d models of plasmas [12] . 

As far as the canonical and the microcanonical ensemble are concerned, a 
comprehensive theory of equivalence and nonequivalence has been worked out (see 
[TBI U\ for a short account and [15] for an extensive treatment). In the following we 
will restrict the discussion to the thermodynamic limit of systems in the microcanonical 
and in the canonical ensemble. Before entering this subject let us recall some basic 
definitions and fix some notation. 

2. 1 . Thermodynamic functions 

We consider Hamiltonian systems with N degrees of freedom, characterized by a 
Hamiltonian function 7i : An — > R of the form 



configuration space and by A N its cotangent bundle, i.e., the phase space. We shall 
usually denote configurations as q = (qi, . . . ,q N ) e T N and phase space points as 
(p,q) = (pi, ...,p N ,q u ...,q N )E A N . 




(1) 



n=\ 



with some potential energy V : T N 



R. We denote by Tn the iV-dimensional 
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The fundamental quantity of the microcanonical ensemble is the Boltzmann entropy 
or microcanonical entropy as a function of the energy (per degree of freedom) e, 

s(e)= lim —In / dpdq5 [H(p, q) - Ne] , (2) 

TV-^oo TV J An 

where 5 denotes the Dirac distribution^ A related quantity is the configurational 
microcanonical entropy as a function of the potential energy (per degree of freedom) 

v, 

s c (v) = lim 4 ln / dqS[V{q)- Nv). (3) 



N^oo N J v 



The fundamental quantity of the canonical ensemble is the canonical free energy 
<p{{3) = - lim ^- In / dpdqe- pn{p ' q \ (4) 



.V 



JV^oo TV J Af 

where (3 = T _1 is the inverse temperature. A related quantity is the configurational 
canonical free energy 

<p (fl = - lim i-ln / dqe-M<). (5) 

Note that the quantity usually called "canonical free energy" in physics textbooks is 
f = if/ (3 and has the dimension of energy (or temperature, since h& — 1). In contrast, 
ip is an adimensional quantity, typically used in the mathematical literature and in 
particular in the recent literature on ensemble (non) equivalence |13[ [14"]. 

By means of Laplace's method for the asymptotic evaluation of integrals [TB] , <p is 
found to be the Legendre-Fenchel transform s* of s, 

<p(P)=s*(P)=M\J3e-s(e)] (6) 

e 

(and analogously for the configurational quantities). 

We may also use the Legendre-Fenchel transform in order to characterize the 
concavity of the entropy A concave function is a real function g : / — ► R defined 
on some interval I for which the inequality 

g(ax + (1 - a)y) ^ ag{x) + (1 - a)g{y) (7) 

with a G [0, 1] holds for all x,y G /. If the inequality is strict, g is called strictly concave. 
The concave envelope of a function s is defined as the smallest concave majorant of s. 
One can show that the concave envelope of a continuous function is given by its twofold 
Legendre-Fenchel transform [T7j 

s**(e)=ip*(e)=M[(3e-tp(P)]. (8) 

A concave function coincides with its concave envelope, whereas a nonconcave function 
does not, 

s(e) is concave <^=^> s(e) = s**(e). (9) 

| We define all thermodynamic functions per degree of freedom, which accounts for the factor 1/TV in 
the definitions. The Boltzmann constant fee is set to unity. 
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This has a remarkable consequence: The free energy tp, by definition, can always be 
obtained from the entropy by means of a Legendre-Fenchel transform. The converse, 
however, is not always true, since the correct entropy s can be computed as the Legendre- 
Fenchel transform of tp only in case s is a concave function. 



2.2. Full equivalence, partial equivalence, and nonequivalence of ensembles 

(Non) equivalence of ensembles can be defined on different levels of description. Here we 
consider ensemble equivalence on the thermodynamic level as discussed in [13], whereas 
a related definition can be given making use of the concept of macrostates p3l H5] . At 
the thermodynamic level one compares the energy e in the microcanonical ensemble 
with the average energy in the canonical ensemble, 

tip) = ^ (n) p , (io) 

where (•}* stands for an average over the canonical probability measure. This leads us 
to the following definitions: 

Full equivalence of ensembles: for any value of e there exists exactly one value of /3 
such that e = (p'(0). 

Partial equivalence of ensembles: for any value of e there exists a value of (3 such 
that e = <p r (P), but the correspondence is not one-to-one. 

Nonequivalence of ensembles: there exist values of e which cannot be obtained as 
canonical average energies for any value of /3. 

Note that this definition of partial equivalence differs slightly from the definition used by 
Ellis, Haven, and Turkington in [15], where only a many-to-one relation between energies 
and temperatures is called partial equivalence, but not the inverse situation of a one-to- 
many relation. From a mathematical point of view, their non-symmetric definition may 
appear reasonable, but from a physical point of view we regard the above definition as 
preferable which treats the microcanonical and the canonical ensemble in a symmetric 
way. 

The above definitions of (non) equivalence can be rephrased in terms of so-called 
supporting lines, and the same can be done for other concepts based on the mathematics 
of Legendre-Fenchel transforms. Supporting lines are a geometric interpretation of the 
Legendre-Fenchel transform (jHJ), and we will use it to some extent implicitly in the 
following. A detailed account can be found in [18] . 



It is useful to translate the properties given in the above definitions of (nonequiv- 
alence into properties of the microcanonical entropy function s. To this aim, we 
distinguish several cases in the following. We restrict ourselves to the case s'(e) > 0. 
First, because it is the physically interesting case of positive temperatures, and second, 
because allowing s'(e) to be negative would require the discussion of further cases 
than those listed below. Rigorous formulations and proofs of several of the following 
statements can be found in [15]. 
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(a) If s is a concave function and its derivative s' takes on all values from the interval 
(0, +00), then equation ([6]) reduces to a standard Legendre transform 

<p(fl = s*((3) = - s(e(/3)), (11) 

where e(/3) is defined implicitly as the solution of the equation 

= /3. (12) 

Geometrically speaking, e{(3) is given as the energy value at which a straight line 
with slope j3 is tangent to the graph of s (see figure [U for an illustration). From 
the above specified properties of the entropy (concavity of s and codomain (0, +00) 
of s'), the one-to-one correspondence of microcanonical and canonical descriptions 
follows and we have full equivalence of ensembles. 



s e 




Figure 1. Sketch of a case of equivalence of ensembles, illustrating the correspondence 
between the microcanonical entropy s(e) and the canonical free energy tp(@). 

(b) If s is not concave, the Legendre- Fenchel transform 

<p(P) = 8*(0)=M\/3e-8(e)] (13) 

£ 

can again be interpreted geometrically: e{(3) is obtained as the energy value at 
which the topmost line with slope f3 tangent to the graph of s touches the graph 
(see figure [2] (top) for an illustration). This implies that, due to the nonconcavity 
of s, there is an interval (si, £2) of energies which do not correspond to any value of 
(3, and therefore the microcanonical and the canonical ensemble are nonequivalent. 
The interval (£1,62) is given by the values of e for which s(e) and its concave 
envelope s**(e) disagree. 

Nonconcave entropy functions may occur in systems with long-range interactions, 
and the nonconcave region of s indicates the presence of a discontinuous phase 
transition in the canonical ensemble. The first solvable model exhibiting nonequiv- 
alence of ensembles was proposed by Hertel and Thirring [19]; other examples are 
the mean-field Blume-Emery-Griffiths model [I] and the mean-field /c-trigonometric 
model for k > 2 [20]. 

(c) If s is concave, but not strictly concave, there exists an interval of energies (£1,62) 
such that s(e) is affine for all e G (£i,£ 2 ), i.e., its graph is a straight line in this 
interval (see figure [2} bottom). Then all the values of the energy in that interval 
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Figure 2. Sketch of nonequivalence (top) due to a nonconcavity of the microcanonical 
entropy s(e), and of partial equivalence (bottom) due to the presence of an affine part 
in the entropy. In both cases the canonical free energy ip((3) is singular and there is a 
discontinuous phase transition in the canonical ensemble. 

correspond to the same value of the inverse temperature P, equal to the slope of 
the affine part of s. Hence, although every value of e corresponds to some value of 
ft, the correspondence is not one-to-one, and we identify this situation as partial 
equivalence of the microcanonical and the canonical ensemble. 
An affine region of s indicates the presence of a discontinuous phase transition 
in the canonical ensemble. In contrast to nonconcave entropy functions, such 
concave, but not strictly concave, entropy functions can occur for systems with 
short-range interactions. In fact, partial equivalence of ensembles is the typical 
situation observed when a discontinuous phase transition takes place in a short- 
range system. 

Other situations of partial equivalence of ensembles occur for concave s when the 
derivative s' does not take on all values in the interval (0, oo). 

(d) Assume that s is a concave function such that s' takes on positive real values, but 
with some interval excluded, 

s'(e) G [0,A]U[/3 2 ,+oo), (14) 

with p 2 > Pi > 0. This happens when the graph of s(e) has a cusp for some value 

e = e c , 

s'{e-) = p 2 > Pi = s'(e+) (15) 

(see figure [3] for an illustration). Then, in the microcanonical ensemble, the inverse 
temperature P(e) = s'(e) jumps discontinuously at e c from j3 2 to Pi, and this 
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behaviour may be termed a microcanonical discontinuous phase transition [81 I2"T|. 
In the canonical ensemble, one obtains a free energy (p(ft) which is affine for 
ft G (fti,ft2): all the values of the inverse temperature in the interval (ftiifh) 
correspond to the same value e c of the energy, and we identify this situation as 
partial equivalence of the microcanonical and the canonical ensemble. 




Figure 3. Sketch of partial equivalence due to a cusp in the microcanonical entropy 
s(e). The canonical free energy ip(/3) has an affine part and there is a discontinuous 
phase transition in the microcanonical ensemble. In the canonical ensemble the 
first derivative of the free energy (the average energy) is continuous, but the second 
derivative has two discontinuities. 

We observe in figure [3] that the singularity of the microcanonical entropy at e c 
corresponds to two singularities in the canonical free energy, located at ft\ and 
fti- For both these inverse temperatures, the first derivative of the free energy is 
continuous, but the second derivative is discontinuous. As a consequence, the heat 
capacity is discontinuous at ft\ and fti, and vanishes for ft 6 (/?i,/3 2 ). Therefore the 
single discontinuous phase transition in the microcanonical ensemble corresponds 
to two continuous phase transitions in the canonical ensemble, and these phase 
transitions take place at the boundaries of the region of partial equivalence. This 
case of partial equivalence is precisely the mirror image of case (c) , with the role of 
the microcanonical and the canonical ensemble interchanged. 
In the physics literature, several examples of systems showing discontinuous 
microcanonical phase transitions are discussed: gravitational systems [T9l |22| 
[23j [23], the mean-field Blume- Emery- Griffiths model |3J, or the mean-field k- 
trigonometric model for k > 2 [20]. However, in all these examples — and in contrast 
to the graph sketched in figure [3] (left) — the phase transition point is inside a "region 
of nonequivalence of ensembles" , i. e., inside a range of energies for which the entropy 
s does not coincide with its concave envelope. As a consequence, for the above cited 
systems the respective Legendre-Fenchel transform of s does not show the behaviour 
plotted in figure [3] (right). 

We are not aware of any model from statistical physics properly showing the 
behaviour sketched in figure [3], and one may argue in favour of the claim that 
such a behaviour is impossible to occur. It would be a worthwhile task to make 
this claim more precise and to justify it rigorously. 
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(e) As a final case, assume that the microcanonical entropy s is a concave function 
and that s' takes on values in some interval , oo) with j3 m { n > 0. In particular, 
this situation may occur for a bounded above Hamiltonian function (see figure H] 
for an illustration of this situation). The upper bound of the Hamiltonian 7i gives 
rise to a maximum value e max of the energy (per particle). Under the above stated 
conditions it follows from the Legendre-Fenchel transform ([6]) that, for all values 
P < Aninj the canonical free energy is given by 

(p(P) = const. + f3e max . (16) 

The energy, again for inverse temperatures (3 G [0,/3 min ), is given by 

e(P) = <f/(l3) = e max . (17) 

Hence, we observe that a single energy value e max corresponds to an interval [0, /3 m i n ) 
of inverse temperatures, and we have partial equivalence of the microcanonical and 
the canonical ensemble. This particular situation of partial equivalence has peculiar 




Figure 4. Sketch of partial equivalence with a concave and regular microcanonical 
entropy. Left: microcanonical entropy s(e). Right: canonical free energy <p(/3). 

implications on the analyticity properties of the thermodynamic functions. The 
microcanonical entropy s may be a smooth function, free of any singularities in the 
interior of its domain. Canonically, in contrast, a continuous phase transition at a 
transition inverse temperature /3 m j n is signalled by a nonanalyticity in the canonical 
free energy. 

The situation considered here is of relevance for example for certain spin systems 
(where the potential energy V is typically bounded above and below), and in 
subsequent sections of this article we present two examples showing such behaviour. 

In the above described cases (a)-(e), the role played by the canonical and the 
microcanonical ensembles is not completely symmetric. Nonequivalence of ensembles 
can only be realized as in case (b) by a nonconcave entropy. The mirror image of 
this situation, with the role of s and ip interchanged, is impossible to occur: since the 
canonical free energy (p is given as the Legendre-Fenchel transform of s, it is necessarily a 
concave function. Similarly, the mirror image of case (e) does not exist as a consequence 
of the properties of the canonical free energy (see section 12. ip . 
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From the above discussion, the intimate relationship between (non)equivalence of 
ensembles and phase transitions becomes clear, and for the cases (b) and (c) this has 
been noted a long time ago and discussed in detail pij [131 E!- Other cases like (d) 
and (e) have received considerably less attention, and it is the purpose of the present 
article to partially fill this gap by providing examples of systems which fit into (e), and 
to discuss the physical implications of this type of partial equivalence. 

3. Two systems showing partial equivalence of ensembles 

The two model systems considered in the following are the mean-field XY model and 
the mean-field spherical model. The interparticle interactions of these systems are of 
mean-field- type, i. e., all degrees of freedom are mutually coupled with the same strength. 
Considering, as is often done when studying the equilibrium statistical mechanics of spin 
systems, the configurational thermodynamic functions and (j5J) of these models, we 
find partial equivalence of ensembles in the sense described above. Adding a standard 
kinetic energy term, i.e., a quadratic form in the momenta, full equivalence is restored. 

3.1. Mean-field XY model 

We start by calculating and discussing the configurational entropy and free energy of 
the mean-field XY model in this section. This model is a mean-field version of the XY 
model of magnetism (see e. g. [25J ) . The potential energy is given by 




where the N degrees of freedom $j £ [0, 2tt) are angular variables, yielding as 
configuration space an iV-torus. In order to attain extensivity of the energy, the 
coupling strength between the degrees of freedom is scaled with N. This model has 
been introduced by Antoni and Ruffo in a version with kinetic energy (which we 
will discuss in section H}, and in the literature it is also referred to as "Hamiltonian 
mean-field model". 

The potential energy Vxy is bounded above and below: ^ Vxy ^ y, so 
that the set of the possible values of the potential energy per degree of freedom is 
compact, v £ [0, |]. The model has been solved in the thermodynamic limit — up to 
a maximization procedure over a single variable to be performed numerically — both in 
the canonical and in the microcanonical ensemble J26J [27], and we will state some of 
these results in the following without derivations. 

In the canonical ensemble, for the configurational free energy <p{(3), the expression 





(18) 



2 



max [$(?/)] 



(19) 



y 



is obtained, with 



log[27r/ (j/)] 




(20) 
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(21) 



where Jo is the modified Bessel function of zeroth order. The maximum of $ is found 
by solving the consistency equation for y, 

y_ = h{y) 

where l\ is the modified Bessel function of first order. As long as (3 < 2, the unique 
solution of (EH) is y = 0, yielding a linear behaviour of <p c {(3): the hallmark of partial 
equivalence. Numerically, for (3 > 2 one finds a /3-dependent value y as a solution of the 
consistency equation (12~TT) . corresponding to the maximum in (1191) . Summarizing, in the 
canonical ensemble we have 



f - log(2vr) 







- $ 



for (3^2, 
for (3 > 2, 



(22) 



and the result of a numerical evaluation of this expression is plotted in figure [5] (left). 
The linear part of the graph corresponds to the region of partial equivalence: for (3 < 2, 
the average energy in the canonical ensemble is constant and equal to |. Differentiating 
ip c , one finds 



(v)M = <p'M 




- l- 



(3* 



for (3 < 2, 
for > 2, 



(23) 



and a plot of this quantity can be found in figure [5] (right). The nonanalytic points 
of the configurational canonical free energy and of the canonical potential energy at 
(3 = 2 signal a continuous phase transition from an ordered phase (low temperatures) 
to a disordered phase (high temperatures). 





Figure 5. Configurational canonical free energy <p c {j3) and canonical potential energy 
{v) c ((3) of the mean- field XY model. The vertical lines mark the critical inverse 
temperature f3 = 2. 



A calculation of the configurational microcanonical entropy s c by means of large 
deviation theory has been reported in [27]. As a consequence of the boundedness of 
Vxyi the configurational microcanonical entropy s c is defined on the compact interval 
[0, 2], and one obtains the expression 

s c (v) = -\{v)VT^ + log [h{X{v))l (24) 
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where X(v) is determined by the consistency condition 



I (\(v))VT^2v = h(X(v)). 
The slope of s c is bounded below, 



mm s' c (e) = s' c (1) = 2 
ee[o,i] 



12 



(25) 



(26) 



(see figure [6] for a plot of the graph of s c ). Hence the configurational microcanonical 
entropy displays the features of partial equivalence as specified in case (e) of section 12.21 
and, as discussed there, s c is a smooth function on its domain [0, |]. 




Figure 6. Configurational microcanonical entropy s c (v) of the mean-field XY model. 



3.2. Mean- field spherical model 

As a second example we discuss a mean-field version of the spherical model, introduced 
by Kac and Berlin [28J as a simplified and exactly solvable version of the Ising model of 
ferromagnetism. The potential of the mean-field spherical model is given by 

1 N 

V svh {a) = - — 2J Wj, (27) 
where the N degrees of freedom <7j 6 1R [i = 1, N) are subject to the constraint 

N 

E °i = N - ( 28 ) 

This constraint restricts the space T^v of allowed configurations o = (cr 1 , . . . , cttv) to an 
(iV-l)-sphere with radius y/N \ The potential V sp h is bounded above and below, and 
the possible values of the potential energy per degree of freedom are v G [-§,0]. This 
model has been solved analytically for finite as well as infinite N [29]. Here we are 
interested in the thermodynamic limit. In this case, the configurational microcanonical 
entropy reads 

8c (v) = ^]og(l+2v), (29) 

and the graph of s c is plotted in figure CD s c is a concave function on [— 1,0], and its 
slope is bounded below by s' c (0) = 1. This model therefore represents another example 
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Figure 7. Configurational microcanonical entropy s c (v) of the mean-field spherical 
model. 



of a system displaying partial equivalence as discussed in case (e) of section 12.21 As 
a consequence, the configurational canonical free energy <p c (P) is expected to be affine 
for /3 ^ 1, and this expectation is confirmed by the calculation of the free energy via 
equation (jSJ), yielding [29] 

for /3 < 1, 

+ log/3) for/3>l, 

[see figure [8] (left) for a plot of the graph of <p c ). The canonical potential energy is 

for (3 < 1, 

(v) Q (p) = { 1-/3 „ _ (31) 



(30) 



2/3 



for /3 > 1 



(see figured] (right) for a plot). The nonanalyticities of ip c (/3) and (v) c (ft) at /3 = 1 
signal a continuous phase transition between a low-temperature ordered (magnetized) 
phase (/3 > 1) and a high-temperature disordered phase (/3 < 1). 




(*>c 







Figure 8. Configurational canonical free energy <p c {P) and canonical potential energy 
(v) c ((3) of the mean-field spherical model. 



As in the case of the mean-field XY model, we observe a remarkable differ- 
ence between the microcanonical and the canonical ensemble: the configurational 
microcanonical entropy s c is a smooth function on its domain, whereas in the canonical 
ensemble a phase transition is signalled by a nonanalytic point of the configurational 
canonical free energy. 
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4. Role of the kinetic energy 

Both systems described above, the mean-field XY model and the mean-field spherical 
model, have a continuous configuration space, and they can therefore be endowed with 
a standard kinetic energy T which is a quadratic form in the momenta pi, 

1 N 

T=\Y,ti- ( 32 ) 
i=i 

In the following we will discuss the implications of adding such a standard kinetic energy 
term to these two models. 

It is common belief that adding a standard kinetic energy term to the Hamiltonian 
function of some model has a trivial effect on its thermodynamic behaviour. This 
is correct only as long as canonical thermodynamic functions are considered: in this 
case, the quadratic form in the momenta results in a Gaussian integral in the canonical 
partition function. This integral can be easily performed, yielding a physically irrelevant 
additive constant in the canonical free energy ip(/3). 

We will show in the following that, contrary to the common belief, a standard 
kinetic energy term may have a non-trivial effect on the thermodynamic behaviour of a 
system when the statistical ensembles are not fully equivalent. For the models discussed 
in this article, adding a kinetic energy has the effect of removing partial equivalence and 
of restoring full equivalence of the canonical and the microcanonical ensemble: in the 
presence of a kinetic energy term, the microcanonical entropy s is no longer defined on 
a compact interval, and its slope has no positive lower bound. Furthermore, and in 
contrast to the configurational entropy s c , a nonanalyticity is present in the interior of 
the domain of s. 



4-1. M ean- field XY model 



Adding a standard kinetic energy term to the potential ([18]) of the mean-field XY model, 
the Hamiltonian function 



1^ 



(33) 



is obtained, where the p^ i = i?j are the momenta conjugate to The dynamical 
properties of this model have been studied in some detail, and peculiar features like 
quasi- stationary states were observed [30J. Calculations of thermodynamic functions of 
this model have been reported for example in [27], where by means of a large deviation 
computation the microcanonical entropy 

jlog« + ac (e-l 

is obtained, where s c is the configurational microcanonical entropy (121 
in (134"|) is attained when 

1 . / u s 



S[E 



sup 

M>0 



(34) 

The supremum 
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and the microcanonical entropy can be written as 

s (e) = ±logu(e) + s c (e-^J, (36) 

where u(e) is given as the solution of ( l35i) . The entropy s is a strictly concave function on 
its domain [0, +oo) and s' takes on all values in [0, +oo). Therefore the microcanonical 
and the canonical ensembles are fully equivalent (see figure [9] (left) for a plot of the graph 
of s). A phase transition is signalled by a nonanalyticity of s(e) at e — § , resulting in 
a kink in s' and in a discontinuity in s" at this value of the energy (see figure [9] (right) 
for a plot of the graph of s'). 




e 



Figure 9. Microcanonical entropy s(e) (left) and its derivative s'(e) (right) for the 
mean-field XY model. 



4-2. Mean- field spherical model 

Adding a kinetic energy term to the potential of the mean-field spherical model, 
the Hamiltonian function 

1 N 

T-C sp h(or,Pa) = 2 Jjftri + "^sph(» (37) 
i=i 

is obtained, where the p ai = &i are the momenta conjugate to <Tj. As a consequence of 
the spherical constraint (125]) . the momenta are subject to the condition 

1 d N N 

=2dt5>* = 5>^- (38) 

i=l i=l 

For the mean-field spherical model, analytic calculations are feasible in both the 
microcanonical and the canonical ensemble, not only in the thermodynamic limit, 
but also for any finite number of degrees of freedom. The computation of the 



microcanonical entropy at finite and infinite is reported in Appendix B In the 
thermodynamic limit, the microcanonical entropy can be written as 

flog [1(1 + 2e)l fore^i, 
s (e) = < , L2 J f (39) 

V ; \|log(2e) for £ >i V ; 

s is a strictly concave function on its domain [— l,oo), and s' takes on all values in 
[0, +oo) (see figure [TUlfor a plot of the graph of s). Hence the microcanonical entropy 
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s of the mean-field spherical model with kinetic energy fulfills the conditions of case 
(a) in section I2.2[ guaranteeing full equivalence of the microcanonical and the canonical 
ensemble. Furthermore, s(e) has a nonanalytic point at £ = |, signalling the occurrence 
of a continuous phase transition at this value of the energy. 



s 




Figure 10. Microcanonical entropy s(e) of the mean-field spherical model with kinetic 
energy term. 

5. Discussion and conclusions 

We have discussed two models, the mean-field XY model and the mean-field spherical 
model, whose configurational microcanonical entropies s c are concave functions defined 
on a compact interval and whose slope s' c has a positive lower bound. As a consequence, 
the microcanonical and the canonical ensemble are partially equivalent. Partial 
equivalence is removed and full equivalence is restored when, instead of configurational 
quantities, the "full" thermodynamic functions s and tp are considered. This behaviour 
has remarkable consequences, which are the main results of the present article and which 
are summarized and discussed in the following. 

At variance with the common belief, the contribution of a standard kinetic energy 
term to the thermodynamic behaviour of a system is not always trivial. Whereas 
in the canonical ensemble the addition of a kinetic energy term results only in a 
physically irrelevant additive constant in the free energy <p(/3) (corresponding to a linear 
contribution to the "usual" free energy /(T)), the effect of such a term can be more 
substantial in the microcanonical ensemble. For the models discussed in the present 
article, we have observed the following effects when adding a kinetic energy term to the 
respective potential energy functions: 

(i) Partial equivalence of the microcanonical and the canonical ensemble is removed 
and full equivalence is restored. 

(ii) The configurational microcanonical entropy s c is a smooth function for both, the 
mean-field XY model and the mean-field spherical model. Adding a standard 
kinetic energy term, a nonanalyticity shows up in the microcanonical entropies s, 
signalling the occurrence of a phase transition in these models. This indicates that 
care is advisable when, in order to simplify a computation, a kinetic energy term 
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is disregarded and only the configurational part of some model is considered: in 
order to identify a phase transition, it is in general not sufficient to consider the 
analyticity properties of the configurational microcanonical entropy, but special 
attention has to be paid to the influence of the boundary points of its domain. 

The conclusions drawn from the two models discussed should extend to all systems with 
bounded above potential energy function V for which the slope of the configurational 
microcanonical entropy s c has a positive lower bound (case (e) of section [2T2T) . 

So for which type of models can this happen? Boundedness of potential energy 
functions is a widespread property among spin models, but only a few examples are 
known for which the slope of the configurational microcanonical entropy has a positive 
lower bound. To the best of our knowledge, all these examples are long-range models, 
but we are not aware of any argument confirming the necessity of long-range interactions 
to obtain such a bound of s' c . This might be an interesting point for future investigations. 
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Appendix A. Mean-field XY model with external magnetic field 

Partial equivalence of the microcanonical and the canonical ensemble can be removed 
not only by adding a kinetic energy term. In this appendix, it is shown that partial 
equivalence may be removed and full equivalence may be restored also by switching on 
an arbitrarily small external magnetic field. 

Adding to the potential ( fTHl) of the mean-field XY model discussed in section 13.11 
a coupling to an external magnetic field h, 



the configurational microcanonical entropy can again be computed by a large deviation 
technique. We refrain from writing down the analytic expression for s c in this case, but 
present a plot of the graph in figure |A~T1 As may be read off from the graph, the slope of 
s c is not anymore bounded, and instead of partial equivalence we find full equivalence of 
ensembles for the mean-field XY model in the presence of a non-zero external magnetic 




(A.l) 



field. 
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Figure Al. Configurational microcanonical entropy s c (v) of the mean-field XY model 
with external magnetic field h = i. Note that the slope of s c is unbounded above and 
below. 



Qualitatively the same happens for the mean-field spherical model: adding a 
coupling to an external magnetic field to the potential ff27|) of the mean-field spherical 
model, we obtain 

^ N N 



2N ^ 3 



i=l 



where the <7j are still subject to the constraint (|28|) . The configurational microcanonical 
entropy for this model is given by 



s c (v) 



In 



1 - (\h\ - Vh 2 - 2v 



(A.3) 



(see [29J for a derivation of this result and for a plot of the graph of s c ). Like in the case 
of the mean-field XY model, the slope of s c has not anymore a positive lower bound, 
and we observe full equivalence of the microcanonical and the canonical ensemble. 



Appendix B. Mean-field spherical model in the microcanonical ensemble 

The fundamental quantity of the microcanonical ensemble is the density of states as a 
function of the energy per degree of freedom e which, for a system with iV degrees of 
freedom characterized by a Hamiltonian function : R, is defined as 



f2Ar(e) oc 



/ dx5[H N {x)-Ne). 
Ja n 



(B.l) 



We restrict ourselves to proportionalities, since multiplicative constants in the density 
of states result only in physically irrelevant additive constants in the microcanonical 
entropy sn = In Qn- For the mean- field spherical model as defined by the Hamiltonian 
function (1371) . the density of states can be written as 

2 



Otv(^) oc / da da 5 



N 



i=l 



8 



N 



i=l 



5 




N 



°i + 2N ? 



i=i 



(B.2) 
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where a = (ci, . . . , <Jn) t and a = ■ ■ ■ , &n) t ^ The first ^-distribution in this 
expression accounts for the spherical constraint (1281) . the second ensures the velocity 
cr to be tangent to the configuration space of the spherical model (constraint ( |38|) ). and 
the third foliates phase space into shells of constant energy Ne. The form of the integral 
in (1B.2I) suggests the use of iV-dimensional spherical coordinates r, (p, . . . and f, (p, . . . 
for the a and the a integrations, respectively. With this choice of coordinates, we have 
J2iLi °1 = r<2 an d J2iLi^i = f 2 - Choosing the polar axis of the position coordinates 
cr along the space diagonal (1, . . . , 1), and the polar axis of the momentum coordinates 
& in the direction of a, we can write Yli=x (T i = V~Nr cos $ and X)i=i a i&i = rr cos , d. 
(Note that the first sum is the scalar product of a and (1, . . . , 1), the second sum is the 
scalar product of a and &.) These orientations of the polar axes render the integrand 
of ( 1B.2j) independent of all angles but the polar angles $ and Performing the trivial 
2(N — 2) other angle integrations, we obtain the proportionality 

p oo roc pir 

tt N (e)oc / drr"" 1 5 (r 2 - N) / dfr^ 1 / d$ sm N ~ 2 $ 
Jo Jo Jo 

dtf sm N ~ 2 i) 5 (rr cos tf) 5 (r 2 cos 2 i? - r 2 + 2Ne) . (B.3) 

The first and second of the ^-distributions allow to easily perform the r- and the •&- 
integrations, yielding 



Sl N (e) oc J d$ sm N ~ 2 ■& dr r N ~ 3 5 (f - ^(cos 2 d + 2efj . (B.4) 



We observe that this integral can be non-vanishing only if cos 2 $ + 2e ^ 0, and, 
performing the r-integration, we write 

n N (e) oc f d$ sin N - 2 (cos 2 # + 2e) {N ~ 3)/2 9 (cos 2 + 2e) , (B.5) 

where 9 denotes the Heaviside step function. With the substitution of variables 
y = cos 2 we obtain the expression 

Q N (e) oc / dyy- 1 ' 2 (1 - y)^ 3 )/ 2 (2e + yf^ 2 Q (2e + y) (B.6) 
>/ o 

for the density of states. Distinguishing the three cases 

for e ^ — ~, 

(1 _ y )(^ 3 )/ 2 ( 2£ + y)^- 3 )/ 2 for -I < £^ 0, 



Q/vfe) oc < 



-2e 
1 



(B.7) 



dyy-V* (1 - y )^- 3 )/ 2 (2e + y) {N - 3)/2 for < e, 



o 



we observe that the (0 < e)-case is, apart from some prefactor, the standard form of 
the integral representation of the Gauss hypergeometric function 2 i*i [21], whereas the 
integral for — | < e ^ can be transformed into such a standard form by a substitution 

§ To ease the notation, we use a instead of p a to denote the momenta conjugate to a. 
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of variables y — > 1 — y(l + 2e). Then, in terms of Gauss hypergeometric functions and 
Gamma-functions T, we obtain 

{0 for e ^ -§, 

F (^) r (f ) C 1 + 2e )"~ 2 2^1 (i iV - 1; 1 + 2s) for -§ < e < 0, (B 
V^T (N — 1) (25)^- 3 )/ 2 a Fi (§, ^, f ; -i) for < e. 

It follows from the properties of these functions that, in the interior of its domain 
[— |, oo), the density of states fijv(e) has precisely one nonanalytic point at e = for all 
system sizes N ^ 2. 

In order to compute the thermodynamic limit of the microcanonical entropy, 

s(e) = lim llnO^e), (B.9) 

N^oo iv 

we perform a substitution of variables y = x 2 in (IB. 61) . obtaining 

Q N (e) oc f l dx [(1 - x 2 ) (2e + x 2 )] {N ' 3)/2 (B.10) 

Jmax{0,-2e} 

for e > — |. The large- iV-limit of this integral can be evaluated by Laplace's method [16J, 
yielding, to leading order, the maximum of its integrand on the domain of integration. 
Hence, omitting a physically irrelevant additive constant, we can write 

s(e)=lim— In max \(l - x 2 ) (2e + x 2 )] {N ~ 3)/2 . (B.ll) 



The maximum in this expression is attained at x = \ \ — e for — ~ < e ^ |, and at 



x = for e > 2, and we obtain 

rin[i(l + 2e)l fore^i, 

«e = i, ^ f (R12) 

[^ln(2e) fore>i, 

for the microcanonical entropy in the thermodynamic limit. The same result can be 
obtained directly from equation ( IB. 81) by making use of the asymptotic expressions for 
Gauss hypergeometric functions 2 -Pi [32]. The singularity which was present at finite N 
for e = has now moved to e = |. This was previously noted and discussed in reference 
[33] where a short account of the solution of the mean-field spherical model with kinetic 
energy was given. 
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